Search for: All records

Abstract Key functions of antibodies, such as viral neutralisation, depend on high-affinity binding. However, viral neutralisation poorly correlates with antigen affinity for reasons that have been unclear. Here, we use a new mechanistic model of bivalent binding to study  >45 patient-isolated IgG1 antibodies interacting with SARS-CoV-2 RBD surfaces. The model provides the standard monovalent affinity/kinetics and new bivalent parameters, including the molecular reach: the maximum antigen separation enabling bivalent binding. We find large variations in these parameters across antibodies, including reach variations (22–46 nm) that exceed the physical antibody size (~15 nm). By using antigens of different physical sizes, we show that these large molecular reaches are the result of both the antibody and antigen sizes. Although viral neutralisation correlates poorly with affinity, a striking correlation is observed with molecular reach. Indeed, the molecular reach explains differences in neutralisation for antibodies binding with the same affinity to the same RBD-epitope. Thus, antibodies within an isotype class binding the same antigen can display differences in molecular reach, substantially modulating their binding and functional properties. 

We develop a convergent reaction-drift-diffusion master equation (CRDDME) to facilitate the study of reaction processes in which spatial transport is influenced by drift due to one-body potential fields within general domain geometries. The generalized CRDDME is obtained through two steps. We first derive an unstructured grid jump process approximation for reversible diffusions, enabling the simulation of drift-diffusion processes where the drift arises due to a conservative field that biases particle motion. Leveraging the Edge-Averaged Finite Element method, our approach preserves detailed balance of drift-diffusion fluxes at equilibrium, and preserves an equilibrium Gibbs-Boltzmann distribution for particles undergoing drift-diffusion on the unstructured mesh. We next formulate a spatially-continuous volume reactivity particle-based reaction-drift-diffusion model for reversible reactions of the form A + B <–> C. A finite volume discretization is used to generate jump process approximations to reaction terms in this model. The discretization is developed to ensure the combined reaction-drift-diffusion jump process approximation is consistent with detailed balance of reaction fluxes holding at equilibrium, along with supporting a discrete version of the continuous equilibrium state. The new CRDDME model represents a continuous-time discrete-space jump process approximation to the underlying volume reactivity model. We demonstrate the convergence and accuracy of the new CRDDME through a number of numerical examples, and illustrate its use on an idealized model for membrane protein receptor dynamics in T cell signaling. 

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst’s broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation. 

In particle-based stochastic reaction–diffusion models, reaction rates and placement kernels are used to decide the probability per time a reaction can occur between reactant particles and to decide where product particles should be placed. When choosing kernels to use in reversible reactions, a key constraint is to ensure that detailed balance of spatial reaction fluxes holds at all points at equilibrium. In this work, we formulate a general partial-integral differential equation model that encompasses several of the commonly used contact reactivity (e.g., Smoluchowski-Collins-Kimball) and volume reactivity (e.g., Doi) particle models. From these equations, we derive a detailed balance condition for the reversible A + B ⇆ C reaction. In bounded domains with no-flux boundary conditions, when choosing unbinding kernels consistent with several commonly used binding kernels, we show that preserving detailed balance of spatial reaction fluxes at all points requires spatially varying unbinding rate functions near the domain boundary. Brownian dynamics simulation algorithms can realize such varying rates through ignoring domain boundaries during unbinding and rejecting unbinding events that result in product particles being placed outside the domain. 

Protein–protein binding domains are critical in signaling networks. Src homology 2 (SH2) domains are binding domains that interact with sequences containing phosphorylated tyrosines. A subset of SH2 domain–containing proteins has tandem domains, which are thought to enhance binding affinity and specificity. However, a trade-off exists between long-lived binding and the ability to rapidly reverse signaling, which is a critical requirement of noise-filtering mechanisms such as kinetic proofreading. Here, we use modeling to show that the unbinding rate of tandem, but not single, SH2 domains can be accelerated by phosphatases. Using surface plasmon resonance, we show that the phosphatase CD45 can accelerate the unbinding rate of zeta chain–associated protein kinase 70 (ZAP70), a tandem SH2 domain–containing kinase, from biphosphorylated peptides from the T cell receptor (TCR). An important functional prediction of accelerated unbinding is that the intracellular ZAP70–TCR half-life in T cells will not be fixed but rather, dependent on the extracellular TCR–antigen half-life, and we show that this is the case in both cell lines and primary T cells. The work highlights that tandem SH2 domains can break the trade-off between signal fidelity (requiring long half-life) and signal reversibility (requiring short half-life), which is a key requirement for T cell antigen discrimination mediated by kinetic proofreading. 

For a chemical signal to propagate across a cell, it must navigate a tortuous environment involving a variety of organelle barriers. In this work we study mathematical models for a basic chemical signal, the arrival times at the nuclear membrane of proteins that are activated at the cell membrane and diffuse throughout the cytosol. Organelle surfaces within human B cells are reconstructed from soft X-ray tomographic images, and modeled as reflecting barriers to the molecules’ diffusion. We show that signal inactivation sharpens signals, reducing variability in the arrival time at the nuclear membrane. Inactivation can also compensate for an observed slowdown in signal propagation induced by the presence of organelle barriers, leading to arrival times at the nuclear membrane that are comparable to models in which the cytosol is treated as an open, empty region. In the limit of strong signal inactivation this is achieved by filtering out molecules that traverse non-geodesic paths.